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Abstract 


The  Biot  equations  are  the  most  promising  approach  for  modeling 
the  propagation  of  acoustic  waves  in  ocean  sediments.  A  method  has 
been  developed  by  Bedford,  Oostley  and  Stem  for  determining  the 
drag  and  virtual  mass  coefficients  in  Biot’s  equations  as  functions  of 
frequency.  The  method  requires  solving  for  the  motion  of  the  fluid  in 
the  pores  when  the  pore  walls  are  subjected  to  a  spatially  uniform, 
oscillatory  motion.  In  this  project,  a  finite  element  algorithm  has 
been  developed  to  determine  the  fluid  motion.  The  accuracy  of  the 
model  has  been  verified  using  a  model  problem  for  which  exact  results 
could  be  obtained.  The  drag  and  virtual  mass  coefficients  have  been 
determined  for  several  two  dimensional  pore  spaces.  It  has  been  de¬ 
termined  that  the  drag  coefficient  is  insensitive  to  the  pore  geometry 
while  the  virtual  mass  coefficient  is  very  sensitive  to  the  pore  geom¬ 
etry.  It  has  also  been  shown  that  the  coefficients  can  be  expressed 
in  nondimensional  forms  which  permit  their  values  to  be  scaled  for 
different  values  of  a  characteristic  linear  dimension  of  the  pore  space. 
This  project  has  established  a  firm  basis  for  proceeding  to  determine 
the  coefficients  for  three  dimensional  pore  spaces. 


Approved  for  public  re?  .- 
Platdbubon  •Unlimited 


j  -L.  ' 

\I 


2J 


87 


Biot  (1956a,  1956b)  developed  a  theory  to  model  the  propagation  of  acous¬ 
tic  waves  in  fluid  saturated  porous  media.  It  has  been  shown  by  numerous 
investigators  that  the  Biot  theory  correctly  predicts  many  of  the  acousti¬ 
cal  properties  of  fluid  saturated  porous  media;  see,  for  example,  Berryman 
(1969),  Hovem  &  Ingram  (1979),  Johnson  &  Plona  (1982),  Stoll  (1977), 
and  Stoll  k  Bryan  (1969).  As  a  result,  this  theory  is  regarded  by  many 
investigators  to  be  the  most  promising  avenue  for  the  study  of  waves  in 
saturated  porous  media. 

However,  an  important  question  remains  to  be  answered.  It  is  known 
that  the  coefficients  in  Biot’s  equations  depend  upon  the  wave  frequency, 
and  they  also  depend  on  the  fabric,  or  microstructure,  of  the  solid  con¬ 
stituent.  Until  these  coefficients  can  be  determined  as  functions  of  fre¬ 
quency,  an  important  obstacle  to  the  accurate  comparison  of  the  Biot  the¬ 
ory  with  acoustic  measurements  in  porous  media  will  remain. 

Recently,  Bedford,  Costley,  &  Stern  (1984)  proposed  a  technique  for 
evaluating  the  drag  and  virtual  mass  coefficients  in  Biot’s  equations  as 
functions  of  frequency.  The  method  requires  solving  for  the  motion  of  the 
fluid  in  the  pores  when  the  pore  walls  are  subjected  to  a  spatially  uniform, 
oscillatory  motion.  To  evaluate  the  method,  Bedford  (1986)  has  applied 
it  to  a  medium  consisting  of  alternating  plane  layers  of  fluid  and  solid. 
For  this  simple  medium,  the  drag  and  virtual  mass  coefficients  could  be 
determined  in  closed  form.  Furthermore,  in  this  case  the  phase  velocity 
and  attenuation  of  plane  waves  could  be  determined  exactly. 

It  was  found  that  when  the  drag  and  virtual  mass  coefficients  were 
determined  using  the  method  of  Bedford,  Costley  &  Stem,  the  velocity 
and  attenuation  of  the  fast  and  slow  waves  predicted  by  the  Biot  theory 
agreed  extremely  closely  with  the  first  two  modes  of  the  exact  solution  over 
a  large  range  of  frequencies.  However,  when  the  virtual  mass  coefficient 
was  assumed  to  be  independent  of  frequency  (which  is  currently  done  by 
most  investigators),  the  Biot  theory  accurately  predicted  the  velocity  and 
attenuation  only  at  very  low  frequency. 

With  these  results  as  motivation,  a  program  of  research  has  been  initi¬ 
ated  to  use  the  finite  element  method  together  with  the  method  of  Bedford, 
Costley  &  Stern  to  determine  the  Biot  drag  and  virtual  mass  coefficients 
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for  realistic  pore  spaces.  In  this  paper,  the  finite  element  formulation  is 
described  and  used  to  determine  the  coefficients  as  functions  of  frequency 
for  various  two  dimensional  pore  spaces. 


The  Method 


The  one  dimensional  forms  of  Biot’s  equations  can  be  written  (Biot  1956a) 


(i -ov.fi.  = 

-  iij)  -  c(u,  -  uf )  +  /,, 


d2u. 


>d2u, 


(1) 


Wm  =  Q-w+R-dx, 


+£(«,  -  iif)  +  c(u,  -  iif)  +  //. 

The  first  equation  is  the  equation  of  motion  of  the  solid  constituent.  The 
term  <f>  is  the  porosity  (the  pore  volume  per  unit  volume  of  the  porous 
medium).  The  terms  p„  and  ua  are  the  density  and  the  displacement  of  the 
solid  material.  The  terms  P,  p,  Q  and  R  are  constitutive  coefficients  which 
depend  on  the  properties  of  the  fluid  and  solid  constituents.  The  terms 
containing  b  and  c  are  forces  exerted  on  the  solid  constituent  by  the  fluid 
constituent  due  to  their  relative  motion.  The  term  containing  b  is  linear  in 
the  relative  velocity,  and  b  is  the  drag  coefficient.  The  term  containing  c 
is  linear  in  the  relative  acceleration.  This  is  called  a  “virtual  mass”  force, 
and  c  is  the  virtual  mass  coefficient.  The  second  equation  is  the  equation 
of  motion  of  the  fluid  constituent.  The  terms  pj  and  Uf  are  the  density  and 
the  displacement  of  the  fluid.  The  terms  f,  and  //  are  external  body  force 
densities. 

Consider  an  imaginary  experiment  in  which  the  solid  constituent  is  given 
a  spatially  uniform,  steady  state  oscillatory  motion 


u,  =  Deiwt, 


(2) 


where  D  is  a  real  constant  and  u  is  the  frequency.  This  could  be  achieved 
(hypothetically)  by  applying  a  suitably  prescribed  external  body  force  den¬ 
sity  ft-  The  resulting  steady  state  motion  of  the  fluid  will  be  of  the  form 

uj  =  Ueiwt,  (3) 


.•j:Uy  Ccxles 
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where  U  is  a  complex  constant.  Substituting  these  two  expressions  for  us 
and  uj  into  the  Biot  equation  of  motion  for  the  fluid  constituent  (1)2  gives 
the  result 

4>PjuU  =  (ib  —  wc)(U  -  D).  (4) 

The  terms  in  the  equation  (1)2  containing  Q  and  R  vanish  because  the 
motion  is  spatially  uniform.  The  complex  equation  (4)  can  be  solved  for 
the  two  real  coefficients  b  and  c,  yielding 

b  =  4>pju  Imaginary!# /(#  -  1)],  (  . 

c=  -0p/Real[#/(#-l)],  1  ' 

where  U  =  U/D.  If  the  term  U  is  known  as  a  function  of  frequency,  Eqs.  (5) 
can  be  solved  for  b  and  c  as  functions  of  frequency.  Thus  determining  6  and 
c  reduces  to  the  solution  of  a  steady  state  boundary  value  problem  in  fluid 
mechanics:  give  the  boundary  of  the  pore  volume  of  a  porous  medium  a 
spatially  uniform,  oscillatory  motion  (2)  and  determine  the  motion  of  the 
fluid  within.  The  volumetric  average  of  the  fluid  displacement  is  identified 
with  the  fluid  displacement  Uf  in  Biot’s  equations.  In  this  way,  U  can  be 
determined  as  a  function  of  frequency. 

Biot  (1956b)  determined  the  drag  coefficient  as  a  function  of  frequency 
by  subjecting  the  pore  fluid  in  a  straight,  cylindrical  pore  to  a  oscillatory 
pressure  gradient.  His  method  was  extended  by  Hovem  &  Ingram  (1979)  to 
determine  the  virtual  mass  coefficient  as  a  function  of  frequency.  Bedford, 
Costley  &  Stern  (1984)  made  numerical  comparisons  of  the  results  of  their 
method  for  the  case  of  a  cylindrical  pore  and  found  that  they  were  identical 
to  those  of  Biot  and  of  Hovem  &  Ingram. 

Can  the  coefficients  determined  as  functions  of  frequency  in  this  way 
be  applied  to  the  problem  of  waves  propagating  in  the  medium?  To  inves¬ 
tigate  this  question,  the  method  has  been  applied  to  a  medium  consisting 
of  alternating  plane  layers  of  fluid  and  solid  (Bedford  1986).  The  material 
properties  and  the  layer  thicknesses  were  chosen  to  correspond  to  water 
saturated  sand.  In  Figs.  1  and  2,  the  phase  velocity  and  attenuation  of 
the  fast  wave  predicted  by  Biot’s  equations  are  compared  to  the  first  mode 
of  the  exact  solution  for  compressional  waves  propagating  parallel  to  the 
layers.  When  the  drag  and  virtual  mass  coefficients  are  determined  using 
the  method  of  Bedford,  Costley  and  Stern,  the  Biot  theory  matches  the 
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exact  solution  very  closely  up  to  a  frequency  of  1  MHz.  (The  curves  are 
indistinguishable.)  However,  when  the  virtual  mass  coefficient  is  assumed 
to  be  constant,  the  Biot  theory  agrees  with  the  exact  solution  only  at  very 
low  frequency. 

Since  the  coefficients  are  determined  by  subjecting  the  solid  constituent 
to  a  spatially  uniform  oscillation,  the  results  will  clearly  be  applicable  to 
propagating  waves  only  when  wavelengths  are  large  in  comparison  to  a 
characteristic  linear  dimension  of  the  pore  space.  However,  this  is  not  an 
important  restriction  since  the  Biot  theory  itself  is  subject  to  the  same 
condition. 

Finite  Element  Formulation 

The  Boundary  Value  Problem 

The  method  of  Bedford,  Costley  &  Stern  requires  that  the  motion  of  the 
viscous,  compressible  fluid  in  the  pore  space  of  a  porous  medium  be  de¬ 
termined  when  the  pore  walls  are  given  a  spatially  uniform,  oscillatory 
motion.  In  this  paper  the  motion  of  the  fluid  is  determined  by  using  the 
finite  element  method.  (See,  for  example,  Becker  et  al.  1981.) 

Under  the  assumption  of  small  displacements,  a  viscous  compressible 
fluid  is  governed  by  the  equation 

P/£i  =  -Vp+  (k+  ^)V(V  -u)  +  rjV2ii.  (6) 

The  vector  u  is  the  displacement  of  the  fluid,  p  is  the  pressure,  k  is  the  bulk 
viscosity,  and  rj  is  the  viscosity.  The  pressure  is  related  to  the  dilatation  of 
the  fluid  by 

p  =  -KV  ■  u,  (7) 

where  K  is  the  bulk  modulus  of  the  fluid. 

The  boundary  condition  imposed  at  the  boundary  of  the  pore  space  is 

u  =  Deiwtcu  (S) 

where  ei  is  a  unit  vector  that  specifies  the  direction  of  the  oscillatory  motion 
of  the  boundary  and  w  is  the  frequency.  The  resulting  steady  state  solutions 


(9) 


for  the  displacement  and  pressure  can  be  written 

u=  ue‘wf , 

P  =  peiujt. 

Substituting  these  solutions  into  Eqs.  (6)  and  (7),  the  steady  state  boundary 
value  problem  reduces  to 

AV2u+£V(V-u)  +  u=0  )  .  „ 

V  ■  u  =  —p/K  )  ma  (10) 

u  =  D  ei  on  dVl 

where  Q  is  the  domain  (the  volume)  in  which  the  solution  is  to  be  obtained, 
dfl  is  the  surface  of  fi,  and  A  and  B  are  defined  by 

A  =  ir)l{p}u), 

B  =  K/(pfU2)  +  *(*  +  \rj)/{pfu). 


The  Variational  Formulation 


Multiplying  Eq.  (10)i  by  a  variation  £u  and  integrating  over  the  domain 
yields  the  variational  expression 

f  [AV2u+ BV(V-u)  +  u]-6xidv  =  0.  (12) 

J  n 

By  integrating  by  parts,  this  equation  can  be  written 

f  [AVu  :  V£u  +  £(V  •  u)(V  •  0u)  —  u  ■  6u]  dv  =  0,  (13) 

J  0 

where  the  notation  Vu :  V^u  denotes  the  product 


Vu:  V«u=SE 


duj  dbuj 

t  *~t  dxi  dxi 


Let  a  functional  ?r(u)  be  defined  by 

ir(u)  =  §•  /  [/IVu  :  Vu  +  B(  V  •  u)(V  •  u)  —  u  •  u]  dv.  (14) 

At  small  frequencies,  the  change  in  density  of  the  fluid  and  consequently 
the  dilatation  V-u  are  very  small.  This  results  in  computational  difficulties 


which  can  be  avoided  by  expressing  the  problem  in  a  form  suitable  for  a 
nearly  incompressible  material  (Herrman  1965).  Equation  (10)2  is  intro¬ 
duced  into  the  functional  (14)  as  a  constraint  by  introducing  a  Lagrange 
multiplier  A.  This  results  in  an  “augmented”  functional  defined  by 

tt*(u,  A ,p)  =  \f  [AVu  :  Vu  —  (B/K)pV  ■  u  —  u  •  u 

Ja  (15) 

+2A(  V  •  u  +  p/K)\  dv. 


By  equating  the  variation  of  this  expression  to  zero,  it  can  be  shown  that 
A  =  —Bp/{2K).  Substituting  this  result  into  the  augmented  functional,  it 
can  be  written 

w*(u,p)  =  ±[  [AVu  :  Vu  -  (2 Bp/K)V  •  u  -  ( B/K2)p 2  -  u  u]  dv.  (16) 


Taking  the  variation  of  this  expression  yields  the  variational  formulation 

Sir*  =  /  [AVu  :  V5u  —  ( B/K)(V  •  u  6p+  pV  ■  6u) 

Jq 

—(BlK2)p6p—  u-  £u]cfv  =  0. 


(17) 


The  Finite  Element  Approximation 

A  finite  element  approximation  is  obtained  by  replacing  the  domain  Q  in 
Eq.  (17)  by  a  discretized  domain  fl/,.  The  objective  is  to  seek  approximate 
solutions  Ufc  and  ph  which  depend  upon  the  mesh  size  of  the  d^cretized 
domain.  The  solution  U/,  is  assumed  to  be  contained  in  a  subspace  of 
H1( ft*),  and  the  solution  ph  is  assumed  to  be  contained  in  a  subspace 
of  where  Hn{ SIk)  denotes  the  set  of  functions  on  f Ij,  whose  nth 

derivatives  are  square  integrable.  Equation  (17)  becomes 

/  [AVuft  :  Vtfu*  -  ( B/I<)(V  •  uh  Sph  +  P/.V  •  Suh) 

J n*  (18) 

-{B/K2)ph  Sph  -  u h  •  «5uh]  dv  =  0. 

In  terms  of  suitable  finite  element  basis  functions  and  0,,the  approxi¬ 
mate  solutions  can  be  written 

pfc  =  E0«A,  (19) 

i  » 
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where  and  Pi  are  the  nodal  point  values  of  «/,  and  ph  respectively. 
Substituting  the  expressions  (19)  into  Eq.  (18)  yields  the  linear  system  of 
equations 

T  t.'  v  1  r  fT  1 

(20) 


'  Kn 

K\2 

■  U  ■ 

.  Kn 

K22 

p 

where 

(I<n )ij=  f  [AWn 

Jtti, 

(K„)h  =  f  [-(B/K) V  ■  4,0,]  dv,  (21 ) 

(Knh  =  [  [-(B/K^GiQjldv 


Equations  (20)  have  been  solved  by  applying  the  essential  boundary  con¬ 
ditions  u  =  Dei  on  the  pore  boundaries  and  periodic  boundary  conditions 
where  the  domains  of  interest  exhibit  periodicity. 


Scaling 


Determining  the  Biot  drag  and  virtual  mass  coefficients  as  functions  of  fre¬ 
quency  by  the  method  described  here  requires  a  substantial  investment  in 
computer  time.  An  important  question  is  whether,  once  the  coefficients 
have  been  determined  for  a  particular  fabric  or  microstructure,  the  re¬ 
sults  can  be  scaled  if  the  characteristic  linear  dimension  of  the  pore  space 
changes.  For  example,  if  the  coefficients  were  determined  for  the  case  of 
spherical  grains  with  a  particular  packing,  can  the  results  be  scaled  to  deter¬ 
mine  the  coefficients  for  spherical  grains  of  the  same  packing  but  a  different 
diameter?  In  this  section  it  is  shown  that  the  answer  to  this  question  is 
yes,  at  least  within  a  particular  range  of  frequencies  and  microdimensions. 

By  substituting  Eq.  (7)  into  Eq.  (6),  the  equation  of  motion  of  the 
viscous  compressible  fluid  in  the  pores  can  be  written  in  terms  of  the  dis¬ 
placement  components  Uk  as 


d2uk 
Pl  dt 2 


=  K 


d2m 

dxidxk 


+  («+  H) 


d2Uj 

dxidxk 


d2uk 

^  dxidii  ’ 


(22) 
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In  terms  of  the  steady  state  solution  Uk  =  Uke,u>t ,  this  is 

2  -  r.  d2Uj  .  .  .  d2u,  .  d2tifc 

-  "  =  Ka^  +  +  +  “"a^r  (23) 

Let  D  be  a  characteristic  linear  dimension  of  the  pore  space,  and  let 
dimensionless  displacements  and  coordinates  be  defined  by 

iik  =  Uk/D,  xk  =  Xk/D.  (24) 

In  terms  of  these  dimensionless  quantities,  Eq.  (23)  can  be  written 


-fi‘={[^]+ihy}^+,'bk 


d2ii 


Uk 


diidxi 


where 


jw  =  — ^2-r,  Re  =  ^,  M-eia?-. 

(K/pi)±  1  «+i>J 


(25) 


(26) 


It  is  seen  that  the  equation  which  governs  the  dimensionless  displacement 
iik  (and  therefore  the  dimensionless  average  displacement  in  the  direction  of 
motion  U )  is  characterized  by  three  dimensionless  groups.  The  dimension¬ 
less  group  M  can  be  recognized  as  a  Mach  number  (the  term  ( K/pj )*  =  c 
is  the  speed  of  sound  in  the  fluid),  and  the  dimensionless  groups  Re  and 
Re'  are  Reynolds  numbers.  This  result  together  with  Eqs.  (5)  for  b  and 
c  leads  to  the  conclusion  that  the  quantities  b/(<j>pfw)  and  cJ{(f)pj )  can  be 
expressed  as  functions  of  M,  Re  and  Re'. 

The  Mach  number  M  is  a  measure  of  the  effect  on  the  solution  of  the 
compressibility  of  the  fluid.  For  frequencies  that  are  sufficiently  low  that 
M 2  1,  the  effect  of  compressibility  can  be  neglected  and  the  quantities 

b/(<t>pfu)  and  c/(<f>pf)  can  be  expressed  as  functions  of  Re  alone.  (The 
analogy  with  gas  dynamics  is  obvious.) 

For  this  reason,  the  principle  results  will  be  presented  in  the  next  sec¬ 
tion  as  plots  of  the  dimensionless  drag  coefficient  b/(<f>p/u)  and  dimen¬ 
sionless  virtual  mass  coefficient  c/{<(>pj )  as  functions  of  the  “dimensionless 
frequency”  Re.  This  makes  the  results  independent  of  the  properties  of  the 
fluid,  and  also  makes  them  independent  of  the  characteristic  linear  dimen¬ 
sion  D  so  long  as  the  restriction  M 2  1  is  satisfied.  It  must  be  emphasized 

that  this  restriction  refers  to  the  frequency  range  within  which  the  results 
can  be  scaled.  The  computations  that  have  been  made  did  account  for  the 
compressibility  of  the  fluid. 


9 


Two  Dimensional  Results 

The  first  objective  was  to  verify  the  finite  element  formulation  and  investi¬ 
gate  the  fineness  of  the  mesh  necessary  to  obtain  accurate  results  for  b  and 
c  as  functions  of  frequency.  In  order  to  do  so,  the  first  case  considered  was 
a  pore  volume  consisting  of  a  plane  layer  of  fluid  parallel  to  the  direction  of 
the  oscillatory  motion.  For  this  case,  an  exact  solution  for  the  fluid  motion, 
and  thus  for  b  and  c,  can  be  obtained  (Bedford  1986). 

An  example  of  the  mesh  used  for  this  case  is  shown  in  Fig.  3.  In  Fig.  4, 
the  computed  values  of  the  dimensionless  displacement  ii  =  u/D  (crosses) 
are  compared  to  the  exact  distribution  across  the  layer  for  a  frequency 
of  100  Hz.  At  this  relatively  low  frequency  the  distribution  is  essentially  a 
Poiseuille  one,  and  the  agreement  of  the  finite  element  solution  with  the  ex¬ 
act  solution  is  excellent.  In  Fig.  5,  the  computed  displacement  distribution 
is  compared  to  the  exact  distribution  for  a  frequency  of  100  kHz.  In  this 
case,  a  relatively  coarse  mesh  was  used  (10  elements  across  the  layer),  and 
the  agreement  with  the  exact  solution  deteriorates  near  the  wall.  In  Fig.  6, 
the  same  case  is  presented  with  a  mesh  of  20  elements  across  the  layer,  and 
the  displacement  distribution  near  the  wall  is  modeled  accurately. 

In  Figs.  7  and  8,  the  values  of  b  and  c  obtained  using  the  computed 
displacement  distributions  for  a  layer  of  fluid  are  compared  to  the  values 
obtained  using  the  exact  distributions.  Computations  are  shown  for  two 
values  of  the  layer  thickness.  These  results  illustrate  the  accuracy  of  the 
finite  element  computations.  In  Figs.  9  and  10,  the  same  results  are  pre¬ 
sented  as  plots  of  b/(u)<f>pj)  and  c/(4>pj)  as  functions  of  the  dimensionless 
frequency  Re.  In  these  two  figures  the  scaling  discussed  in  the  previous 
section  is  illustrated.  The  results  for  three  values  of  the  layer  thickness  fall 
on  a  single  curve. 

An  important  result  that  was  observed  is  that  the  valves  ofb  and  c  were 
approximately  constant  (independent  of  frequency)  up  to  a  value  of  Re  of 
approximately  10. 

Another  case  that  was  considered  was  chosen  to  be  a  two  dimensional 
approximation  of  the  pore  space  in  a  porous  medium  consisting  of  spherical 
grains.  The  geometry  was  a  periodic  array  of  cylinders  with  porosity  <f>  = 
0.365.  The  periodic  “microelement”  and  the  finite  element  mesh  used  for 
this  case  are  shown  in  Fig.  11.  An  important  question  that  was  investigated 
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in  this  two  dimensional  study  was  the  sensitivity  of  the  computed  values  of 
b  and  c  to  the  shape  of  the  “grains”  of  the  porous  medium.  For  this  purpose 
computations  were  also  carried  out  for  three  other  geometries,  each  having 
the  same  porosity  (0.365). 

•  A  periodic  array  of  square  rods.  The  periodic  element  and  mesh  are 
shown  in  Fig.  12.  The  motion  was  horizontal;  thus  the  rods  were 
“diamond  shaped.” 

•  A  periodic  array  of  irregularly  shaped  rods  with  elongation  oriented 
horizontally.  The  periodic  element  and  mesh  are  shown  in  Fig.  13. 

•  The  periodic  array  of  irregularly  shaped  elements  with  elongation 
oriented  vertically. 

The  results  for  b/{u<j>pj)  and  c/{4>pj)  as  functions  of  dimensionless  fre¬ 
quency  Re  for  the  various  shapes  are  shown  in  Figs.  14  and  15.  The  charac¬ 
teristic  dimension  D  was  chosen  to  be  the  minimum  gap  width  between  the 
rods.  The  value  of  D  used  in  the  case  of  the  array  of  square  rods  was  0.0451 
mm.  The  value  of  D  used  for  the  other  shapes  was  0.0224  mm.  From  the 
results  it  is  seen  that  the  drag  coefficient  is  remarkably  insensitive  to  the 
geometry  of  the  pore  space.  (It  is  important  to  recall  that  the  porosity  was 
the  same  in  each  geometry.  In  each  geometry  except  the  square  rods,  the 
minimum  gap  between  the  rods  was  the  same.)  By  contrast,  the  virtual 
mass  coefficient  is  seen  to  be  very  sensitive  to  the  geometry.  However,  for 
both  coefficients,  the  qualitative  dependence  on  the  frequency  of  the  various 
geometries  is  very  similar. 


Conclusions 

It  has  been  found  feasible  to  use  the  finite  element  method  together  with 
the  method  of  Bedford,  Costley  h  Stern  to  obtain  accurate  values  of  the 
Biot  drag  and  virtual  mass  coefficients  as  functions  of  frequency.  Within 
a  particular  range  of  frequency,  the  results  can  be  scaled  to  obtain  the 
coefficients  as  functions  of  frequency  for  other  values  of  a  characteristic 
dimension  of  the  pore  geometry.  The  computed  values  of  the  coefficients 
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are  approximately  constant  (independent  of  frequency)  up  to  a  value  of 
dimensionless  frequency  Re  of  approximately  10. 

When  computations  were  carried  out  for  various  pore  geometries  hold¬ 
ing  the  porosity  fixed,  the  drag  coefficient  was  found  to  be  very  insensitive 
to  changes  in  geometry  while  the  virtual  mass  coefficient  was  found  to  be 
quite  sensitive  to  the  geometry.  However,  the  qualitative  dependence  of 
both  coefficients  on  frequency  was  quite  similar  for  the  various  geometries, 
suggesting  that  it  may  be  possible  to  use  empirical  “shape  factors”  to  ac¬ 
count  for  pore  geometry. 

The  ultimate  objective  of  this  research  is  to  determine  the  Biot  drag  and 
virtual  mass  coefficients  as  functions  of  frequency  for  three  dimensional  pore 
spaces.  The  results  reported  in  this  paper  suggest  that  it  will  be  feasible 
to  do  so. 
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Figure  Captions 

Fig.  1.  Phase  velocity  of  the  Biot  fast  wave  and  the  first  mode  of  the  exact 
solution. 

Fig.  2.  Attenuation  of  the  Biot  fast  wave  and  the  first  mode  of  the  exact 
solution. 

Fig.  3.  A  finite  element  mesh  for  a  plane  layer  of  fluid. 

Fig.  4.  Computed  and  exact  displacement  distribution  at  100  Hz. 

Fig.  5.  Computed  and  exact  displacement  distribution  at  100  kHz  with  10 
elements  across  the  layer. 

Fig.  6.  Computed  and  exact  displacement  distribution  at  100  kHz  with  20 
elements  across  the  layer. 

Fig.  7.  Computed  and  exact  values  of  b  as  functions  of  frequency. 

Fig.  8.  Computed  and  exact  values  of  c  as  functions  of  frequency. 

Fig.  9.  Computed  and  exact  values  of  bfu<f>pj  as  functions  of  Re. 

Fig.  10.  Computed  and  exact  values  of  c/4>pj  as  functions  of  Re. 

Fig.  11.  Periodic  element  and  mesh  for  an  array  of  cylinders. 

Fig.  12.  Periodic  element  and  mesh  for  an  array  of  square  rods. 

Fig.  13.  Periodic  element  and  mesh  for  an  array  of  irregularly  shaped  rods. 

Fig.  14.  The  dimensionless  drag  coefficient  as  a  function  of  dimensionless 
frequency  for  various  pore  geometries. 

Fig.  15.  The  dimensionless  virtual  mass  coefficient  as  a  function  of  dimen¬ 
sionless  frequency  for  various  pore  geometries. 
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NORMALIZED  PHASE  VELOCITY 


Figure  1:  Phase  velocity  of  the  Biot  fast  wave  and  the  first  mode  of  the 
exact  solution 
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Figure  2:  Attenuation  of  the  Biot  fast  wave  and  the  first  mode  of  the  exact 
solution 
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Figure  4:  Computed  and  exact  displacement  distribution  at  100  Hz. 
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Figure  5:  Computed  and  exact  displacement  distribution  at  100  kHz  with 
10  dements  across  the  layer. 


Figure  6:  Computed  and  exact  displacement  distribution  at  100  kHz  with 
20  elements  across  the  layer. 
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Figure  15:  The  dimensionless  virtual  mass  coefficient  as  a  function  of  di 
mensionless  frequency  for  various  pore  geometries. 
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